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Abstract. - Impurity substitution of Zn in La-214 and (Y,Ca)-123 high-T c superconductors 
suppresses T c but does not affect appreciably the onset of the pseudogap phase in the underdoped 
region nor optimal doping or the position of the inferred quantum critical point. Based on a 
1/N expansion of the t — J model we explain these findings as well as the similar dependence 
on a magnetic field in terms of a quantum critical point scenario where a flux phase causes the 
pseudogap. 



The quantum critical point scenario Jl], [|, || represents a popular frame to discuss the phase 
diagram of high-T c oxides. Suppressing superconductivity by strong magnet fields it has been 
found experimentally that there exists a critical hole doping S^ CP at zero temperature which 
separates a metallic state at larger dopings from an insulating state at lower dopings. Strong 
fluctuations of the order parameter related to the insulating phase are thought to suppress the 
density of states for 8 < leading to the pseudogap features in the underdoped region 

and to be instrumental for superconductivity around 8 C and, at higher temperatures, for the 
anomalous properties of the normal state in these systems. 

The microscopic nature of the order parameter of the insulating phase and its fluctuations 
are presently not clear. One obvious choice is antiferromagnetismjjj which occurs at T = 
as long-range ordered phase at zero and small dopings. The corresponding zero temperature 
critical point, however, corresponds to a much smaller doping value than the observed one, 
SQ CP ~ 0.17. A reasonable large 8^ CP has been obtained in Ref. || for a scenario with an 
incommensurate charge density wave (ICDW). In this approach the pseudogap features are 
not directly related to the ICDW order parameter but rather connected to strong d-wave 
superconducting fluctuations sustained by ICDW precursors. Related approaches include 
preformed Cooper pairs where phase coherence is achieved below T c fm or RVB spinon pairing^) 
and the 7r-flux phase||, where charge coherence is obtained by Bose condensation of holons. A 
different proposal has been made in Ref. M . Based on a 1 /N expansion for the t — J model the 
quantum critical point was identified with a transition from the normal to a ci-wave flux state 
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occurring near the observed S^ CP for realistic parameters. In this approach optimal doping 
is determined by the onset of the flux phase and the phase diagram in the underdoped region 
is characterized by the competition between the flux and superconducting order parameters 
both having <i-wave symmetry. 

Recently, several experimental results have been published which may be able to confirm or 
to rule out some of the above approaches. Measurements of NMR spin lattice relaxation rates 
in the presence of magnetic fields up to ~ 15 Tesla did not yield appreciable changes for the 
onset temperature T* of the pseudogap phase whereas the superconducting T c was reduced by 
about 8 K |J. A strong suppression of T c and, at the same time, no change in the pseudogap 
was previously reported in Zn-doped YBa 2 Cu 3 06+x@. High resolution photoemission[]lTf , 
electronic Raman spectroscopy |]l2| , NMR[[l3| and heat capacity data|l4|, [ll| show that T* does 
not merge with T c in the overdoped regime, but vanishes near optimal doping. These findings 
indicate that the pseudogap and the superconductivity are different phenomena and not related 
to the same order parameter. Furthermore, it has been found experimentally [fl6| that the 
lowering of the T c curves in Zn doped Yo.sCao.2Ba2Cu307 ((Y,Ca)-123) and L^-zSr^CuC^ 
(La-214), is concentrated around optimal doping and that the optimal doping itself is not 
shifted. It is the purpose of this Letter to investigate the influence of impurity scattering and 
magnetic fields on the phase diagram calculated in Ref. || and to compare the results with 
the above experimental findings. 

We consider a t-J-V model with N degrees of freedom per lattice site on a square lattice. 
Its Hamiltonian can be written in terms of Hubbard's Y-operators as H 



H 



•jv 2s x - x i + AN 2^ X * X 3 

(ij) (ij) 

p = l...N p,q=l...N 

■4^ 2s x i x j + 2s w { ■ { } 

(ij) a 

p,q—l...N p,g = l...N 



The internal labels p,q... consist of a spin label distinguishing spin up and spin down states and 
a flavor label counting N/2 identical copies of the original orbital, (ij) denotes nearest- neighbor- 
sites. The first three terms represent the t-J Hamiltonian, the last term a screened Coulomb 
interaction appropriate for two dimensions and taken from Ref. Jl7|| . In the following we 
express all energies in units of t. The strength of the Coulomb interaction will be characterized 
by its value between nearest neighbor sites V n . n .. In the limit of large N, the interactions 
become purely instantaneous and H can be diagonalized analytically. In the absence of 
impurity scattering, the coexistence state of superconductivity and a staggered (w, ir) flux 
phase can be obtained from a Nambu representation with 4 states yielding four electronic 
bands with dispersion B 

±£;±(k)-±VK(k)±A] 2 + A(k)2, (2) 



where 



e(k) = ^(k)2+0(k)2. (3) 

Here the momenta k are restricted to the new Brillouin zone which is one half of the original 
one. jj, is a renormalized chemical potential, </>(k) the flux order parameter, A(k) the super- 
conducting gap, and e(k) the one-particle energies in the normal state. Both order parameters 
have d-wave symmetry: 0(k) = <^>[cos(fc x ) — cos(k y )], A(k) = A[cos(k x ) — cos(fcy)]. They are 
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Fig. 1. - Phase diagram of the t-J-V model (see text) taking into account the normal (N), supercon- 
ducting (SC), and the flux phase (FL) as well as the coexisting state of superconductivity and flux 
phase (FL+SC). The normal state-flux phase transition ends in a T — quantum critical point at 
fiQ CP ~ 0.115. The dashed line represents the instability towards the flux phase in the superconducting 
state. 



determined by the self-consistent set of equations: 

<*( k ) = o4-£ J ( k + p)^(p) 



(4) 



A(k) 



t J ( k + p) - ( k + p)]^ A (p) - 



(5) 



where 



4^ tanh 
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2T 
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A(k) 
2E + (k) 
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A(k) 
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tanh 
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2T 



(7) 



The resulting phase diagram, calculated using J = 0.3 and V n , n . = 0.5J, is shown in Fig. 
0. Disregarding superconductivity, the second-order normal state-flux phase transition line 
ends in a quantum critical point, denoted by the black dot, at S^ CP ~ 0.115. We find the 
maximum of T c at essentially the same doping because of the strong competition between flux 
and superconducting phase and also that the flux phase instability is only slightly shifted 
by superconductivity (dashed line). 

Now we are going to investigate how the phase diagram in Fig. [l] is affected by impurity 
scattering. In the simplest approximation, the effects of impurities in the normal state can be 
taken into account by introducing a renormalized frequency fig] 



iuj n + iT- 



(8) 
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where T is a scattering rate, here used as a free parameter proportional to the impurity 
concentration. Throughout the flux phase the self-energy due to impurity scattering is still 
diagonal in the 4x4 Nambu representation because the flux order parameter does not couple to 
the impurities. The constant T in Eq. (||) could be improved in a deeper analysis considering 
effects due to the proximity of the Van Hove singularity E^j. However, the interesting doping 
region for superconductivity is in our model not at all correlated with the Van Hove singularity. 
As a matter of fact, the chemical potential for 8 ~ 0.1 is quite far away from the Van Hove 
singularity, which is at S = in our model. In this situation the band can be assumed in a good 
approximation to be structureless and T as a constant, even in the case of strong potential 
scattering jllj. Experimentally it is known that non-magnetic impurities such as Zn and Al 
induce local moments on neighboring Cu sitesjl^, ^o). Both in a d-wave flux phase and a 
d-wave superconductor random local magnetic moments lead only to renormalizations of the 
frequency. As a result they contribute additively to T in Eq. (||) and thus can be accounted for 
by a proper choice for T. It also has been argued|2(J that strong potential scattering near the 
unitary limit is much more important for the reduction of T c than the scattering from induced 
magnetic moments. 

The occupation number /r of an electronic state with energy e in the presence of impurities 
reads 

* = ' v . 



lUJ n 
n 

i i _ 

= 2 + 2^ Im 



,1 .e + *T\ , (\ .e-iT 



(9) 



where tp denotes the digamma function. In the limit of zero impurity concentration T — > 
and /r(e/T) reduces to the usual Fermi function f(x) = l/(e x + 1). In a similar way, we also 
define a function tanhr by 



tanhr (JL) = \ 



fr ( -j^j - fr 



(10) 



For the determination of the superconducting critical temperature T c , the self-consistent 
set of gap equations can be linearized with respect to A(k). The resulting equations are again 
given by Eqs. (M-M) if the function tanh is everywhere replaced by the function tanhr, defined 



in Eq. (10). The solid lines in Fig. g show numerical results for T c as a function of doping 
6 for different scattering rates T, using J = 0.3 and V n . n . = 0.5J. These curves illustrate the 
suppression of T c with increasing scattering rates T — 0, 2 ■ 10~ 3 , 4 • 10 -3 , and 6 • 10 -3 . The 
corresponding changes in T*, determining the phase boundary between the normal state and 
the flux state, are depicted in Fig. ^ by the grey region. The chosen values for V correspond 
roughly to T ~ 1.0T C at optimal doping, and to T ~ 1.5T C in the strongly underdoped region, 
interpolating between the weak- and the strong-coupling regimes. One important result of 
Fig. |^ is that the flux phase boundary S FL (T) is only slightly shifted by impurities, in spite 
of the strong suppression of the superconducting critical temperature. In particular, the zero 
temperature limit of 5 FL , S FL (0), is almost completely independent of the impurity scattering 
rate. Since in our approach the maximum of T c as a function of doping is essentially determined 
by 8 FL (0) this means that the T C (S) curves shrink to 8 FL (0) with increasing scattering rate 
which is a characteristic feature of Fig. ^. Interpreting Fig. @ in terms of a quantum critical 
point scenario means that the corresponding critical doping 5^ CP is given by 8 FL (0) and that 
§Q CP is almost completely independent of the impurity scattering rate. The curves in Fig. 
are in excellent agreement with the corresponding experimental curves in Zn doped (Y,Ca)-123 
and La-214, given in Fig. 2 of Ref. @. 
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Fig. 2. - Solid lines: Suppression of T c by impurity scattering, calculated with the scattering rates F = 
(squares), 2-10~ 3 (circles), 4-10 -3 (diamonds), and 6 -10 (triangles). Grey region: Corresponding 
variation of the transition temperature T* to the flux state. 



We can gain further insight into our results by the following analysis. We first notice that, 
at least for the above values of T, impurity scattering effects lead to an additional smearing of 
the occupation number in Eq. (^) which can be simulated in a very good approximation by 
an effective temperature: 



Hi) 



As a consequence, the flux phase instability in the presence of impurities is roughly determined 
by S F L (T) ~ 5 FL {T + T). But 8 FL is very weakly dependent on T until T ~ ft, so that the 
T = quantum critical point is not expected to be shifted as long as V < ft holds. 

In order to understand better why S FL (T) is almost independent of T for T < ft, we consider 
the flux phase susceptibility x(T,F, S). The instability line S FL (T) is determined in general 
by the equation: 1 = Jx(T, T, 5), where the dependence of the susceptibility on temperature 
and impurities is given by the factor || 

F v (e) = M(e-fi)/T]-M(-e-ft)/T} _ 
e 

A sketch of the numerator and denominator of Fj? is given in Fig. ^. Continuous and 
discontinuous solid lines represent the numerator for r and T = 0, respectively, the 
dashed line the denominator. Due to the weak variation of the denominator around ft and 
— ft, it is clear that this factor is only slightly affected by a possible smearing due to finite 
temperatures or impurity concentrations, as long as T + T < ft holds. Things arc different 
in the case of the superconducting susceptibility. Here the divergence of the denominator 
coincides with the jump in the numerator, so that even a small smearing leads to a strong 
change in \ and a large suppression of T c . 

We would like to mention that a charge-density-wave susceptibility would contain a similar 



factor as in Eq. (12), so that also in this case impurities will not affect substantially the 
function Fp. However, the CDW order parameter has the symmetry of the underlying lattice, 
i.e., s-wave symmetry. Impurities couple in this case directly to the order parameter and 
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Fig. 3. - Numerator (solid lines) and denominator (dashed line) of Fr(e), Eq. (fl^), in the case of the 
flux instability of the normal state (upper panel). The discontinuous and continuous solid lines refer 
to zero and nonzero effective temperature, respectively. For comparison, the lower panel shows the 
corresponding quantities in the case of a superconducting instability. 



the self-energy due to impurity scattering also acquires non-diagonal elements in addition to 
the diagonal ones described by Eq. (||). As a result, one expects that the charge-density 
wave state is sensitive to impurities and the corresponding quantum critical point and optimal 
doping would be shifted by the impurities, in disagreement with Ref. [fl6| . 

We have considered throughout our analysis a (tt, 7r) flux phase and its competition with 
superconductivity. In Ref. pil it was shown that there is a continuous transition line in the 
T — 5 plane from a commensurate to an incommensurate flux state at low temperatures (the 
term "commensurate" is here used with respect to the lattice periodicity and not, as in Ref. 
[ p2| , with respect to the electronic filling). Taking the incommensurability into account the 
largest onset of the flux phase as a function of doping occurs now at T — with a critical 
doping S^ CP ~ 0.135. However, the boundaries in the phase diagram and, in particular, the 
competition between flux and superconducting phases are not much changed by allowing for 
the incommensurability of the flux phase. Disregarding superconductivity we also have studied 
the influence of impurities on the boundary between an incommensurate flux and the normal 
phase. The resulting width in T* for the scattering rates used in Fig. 2 is very similar to 
that shown in this Figure for the commensurate case. In particular, the change of the critical 
doping at T = was for all T's smaller than 0.01 showing that our Fig. 2 calculated for the 
commensurate case is also valid in the incommensurate case in a very good approximation. 
Also the simplified arguments based on Fig. 3 for the robustness of the flux state in contrast to 
the superconducting state with respect to impurities still apply. The finiteness of the chemical 
potential p, reflects the fact that one-particle states which are also not exactly degenerate 
in energy are involved in forming the flux state. This means that finite difference energies 
|e(k) — e(k + Q)| (Q is the wave vector of the flux phase) associated with a large phase space 
are important which can be characterized by a typical energy 2/2. This explains why both the 
commensurate and the incommensurate flux phase behave in a very similar way with respect 
to impurities. 

Our proposed scenario of a flux quantum critical point is also consistent with the NMR 
measurements of Ref. 0. In that paper, a magnetic field H = 14.8 T was shown to yield a net 
reduction of the superconducting critical temperature of AT C = 7.8 K but no corresponding 
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decrease of the pseudogap temperature T* within the experimental uncertainty of 2 %. In our 
theory the pseudogap and the superconductivity arise from two different mechanisms, so that 
a significant reduction of T c is possible in the absence of a corresponding reduction of T* . The 
predominant effect of a magnetic field on the superconducting phase is a reduction of T c in 
order to balance the free magnetic energy related to the Meissner effect. The decrease of T c 
is linear in H for H <C H c {H c being the critical magnetic field), so that the reduction in T c 
is quite effective even for small magnetic fields. On the other hand, the effect of a magnetic 
field on the flux phase is mainly due to the Zeeman splitting AE = gfisH, where g and fiB 
are the g factor and the Bohr magneton, respectively. This energy is about 20 K for H = 14.8 
TQ, and thus much smaller than the width of the electronic band. If we generalize our order 
parameter cj> in the presence of a magnetic field via <f> —* (/> = [<f>i(fi — AE) + <fii(jl + AE)]/2, 
we obtain an effective susceptibility given by \ — [XT (A — + \i (A + AE)]/2. The Zeeman 
splitting AE = 20 K is much smaller than the energy scale set by the bandwidth W « 0.5 
cV, so that we expect a negligible effect on the flux instability. Moreover the shift of T c will 
be only of order H 2 in the magnetic field. 

We have checked the above arguments by calculating explicitely the change in the instability 
from the normal to the flux state in the presence of a Zeeman splitting AE = 20 K, correspond- 
ing to AE — 6 • 10~ 3 t with t — 0.3 eV. We find a shift in doping of the T = quantum critical 
point of about 4 %, which is quite small. This shift disappears with increasing temperature 
because of thermal smearing which makes the Zeeman splitting ineffective. 

In conclusion, we have analyzed the effects of impurity scattering on the flux phase and 
on its interplay with superconductivity within the framework of a quantum critical point 
scenario. We have found that the transition between the flux and the normal phase is essentially 
unaffected by impurities and magnetic fields, in very good agreement with the experimental 
data. This is especially true at zero temperature where we identify the quantum critical 
point with the transition between the normal and the (incommensurable) flux state. We also 
pointed out that a charge density wave as the origin of the pseudogap phase would directly 
couple to impurities in contrast to the flux order parameter and thus be much more sensitive 
to impurities. 

* Present address: Dipartimento di Fisica, Universita di Roma I "La Sapienza" , P.le Aldo 
Moro 2, 00184 Roma, Italy. 
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